Canonical quantization of constrained theories on discrete space-time lattices 



Cayetano Di Bartolo 1 , Rodolfo Gambini 2 , and Jorge Pullin 3 

1. Departamento de Fisica, Universidad Simon Bolivar, 
Aptdo. 89000, Caracas 1080- A, Venezuela. 
2. Instituto de Fisica, Facultad de Ciencias, Universidad de la Republica 
Igud esq. Mataojo, CP 11400 Montevideo, Uruguay 
3. Department of Physics and Astronomy, Louisiana State University, 
202 Nicholson Hall, Baton Rouge, LA 70803-4001 
(May 29th 2002) 

We discuss the canonical quantization of systems formulated on discrete space-times. We start 
by analyzing the quantization of simple mechanical systems with discrete time. The quantization 
becomes challenging when the systems have anholonomic constraints. We propose a new canonical 
£SJ ' formulation and quantization for such systems in terms of discrete canonical transformations. This 

allows to construct, for the first time, a canonical formulation for general constrained mechanical 
systems with discrete time. We extend the analysis to gauge field theories on the lattice. We 
consider a complete canonical formulation, starting from a discrete action, for lattice Yang-Mills 
theory discretized in space and Maxwell theory discretized in space and time. After completing the 
treatment, the results can be shown to coincide with the results of the traditional transfer matrix 
method. We then apply the method to BF theory, yielding the first lattice treatment for such a 
theory ever. The framework presented deals directly with the Lorentzian signature without requiring 
an Euclidean rotation. The whole discussion is framed in such a way as to provide a formalism that 
would allow a consistent, well defined, canonical formulation and quantization of discrete general 
relativity, which we will discuss in a forthcoming paper. 



> 



O 



I. INTRODUCTION 



Discretizations are used in general relativity mainly in three contexts: a) numerical relativity, b) canonical quantum 
gravity, c) path integral quantum gravity. In the first context the aim of discretizing the equations is to solve them on 
a computer numerically. In the second and third contexts the primary goal is to make finite the infinite dimensional 
computations implied by the field theory nature of general relativity (and eventually to solve the theory numerically 



on a computer). It has long been known in numerical relativity that the set of algebraic equations that results of 
discretizing the Einstein equations (for instance, with externally prescribed lapse and shift) is inconsistent in the 
5h . following sense. The discretized evolution equations produce solutions in the future that do not satisfy the discretized 
^P" constraints, even if one started from initial data that satisfied them. The general attitude towards this problem in 
^ , numerical relativity has been that the aim is to provide an approximate description of the continuum theory and 
therefore the inconsistency is accepted as part of the errors of the approximation scheme. This was clearly discussed 
by Choptuik @, who showed that the constraints were preserved at the same order of discretization as the evolution 
■ equations. 

In canonical and path integral quantum gravity, the potential inconsistency of the discretized theories has not 
received a lot of emphasis. It was noted by Bander M in general and later by Friedman and Jack |3) in the context 
of Regge Calculus that it was not straightforward to find consistent Hamiltonian formulations of general relativity 
discretized in space and with continuum time. It has generally been recognized that a problem has to arise in lattice 
treatments of general relativity since, unlike the case of Yang-Mills theory, the introduction of a lattice breaks the 
"gauge symmetry" of the theory (invariance under diffeomorphisms) . In spite of these observations, the remarkable 
formulation of a well defined Hamiltonian constraint for quantum general relativity due to Thiemann j| starts with 
a discretization of the theory. Thiemann promotes the Hamiltonian constraint to a quantum operator acting on 
the space of solutions of the diffcomorphism constraint. This somewhat obscures the issue of the implementation of 
the constraint algebra, since on such space the algebra of constraints is Abelian. It has been noted, however, that 
it is likely that problems are implicit in the quantization with respect to the closure of the constraint algebra j^] 
"off-shell" (on non-diffeomorphism invariant states). Proposals to construct quantizations of the Hamiltonian that 
produce non- Abelian operators have failed to yield quantities that reproduce the proper constraint algebra 

In canonical formulations, the role of the constraints is to generate infinitesimal symmetries. In a discrete theory, the 
best we could hope for is that the constraints generate some finite version of the symmetry in question. Obviously, such 
an action cannot be achieved with a discretization of the continuum constraints, since the latter generate infinitesimal 
symmetries. The continuum constraints therefore naturally structure themselves into algebras. The discrete generators 
of symmetries that one could have in a discrete theory will not structure themselves in an algebra. Therefore there 
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is no sense in which one will have in the discrete theory an algebra of constraints that "in the limit" reproduces the 
continuum algebra. 

An additional element is that discrete symmetries are not characterized by free parameters, as infinitesimal sym- 
metries are. A discrete transformation can only admit a set of fixed values for the parameters of the transformation. 
Therefore in dealing with finite symmetries the finite quantities that are remnants of the Lagrange multipliers of the 
continuum theory become fixed. This inevitably leads us to systems of constraints of second class in the canonical 
treatment. 

Finally, one has a choice as to the role of time. One could try to discretize "space" and leave "time" continuous. 
This is highly unnatural in general relativity and might further complicate the issues of the constraint algebra since 
the latter mixes symmetries in space and time. In this paper we will choose to pick a discrete time in the theories we 
consider. This is more symmetric and allows for a better connection with "spin-foam" approaches in which time is 
discrete. 

Having a discrete canonical formulation for general relativity is, in our opinion, fundamental at the time of under- 
standing current "spin-foam" approaches to the path integral of quantum general relativity. It is well known that 
path integrals for gauge field theories should be evaluated with care in order to avoid integrating over the orbits of the 
gauge groups. This is at the core of the Fadeev-Popov technique and can only be fully understood if one understands 
the degrees of freedom of the theory at a Hamiltonian level M . To achieve this for discrete gravity is one of the goals 
we are pursuing. 

An additional point to consider is that most lattice treatments of gauge theories require an Euclidean rotation in 
order to be well defined. This is expected to be problematic in the case of general relativity. The formalism we will 
present in this paper does not require an Euclidean rotation and operates directly in terms of the Lorentzian signature 
theory. The lack of a Lorentzian formulation on the lattice has been recognized as one of the major unsolved issues 
in lattice formulations of gravity || 

The inconsistencies that arise in discretizations are present in spin foam approaches directly. When one discretizes 
the action in order to compute the path integral one is left with a discrete action whose equations of motion are 
seldom analyzed and in general yield an inconsistent theory. If one therefore computes the integral for the discrete 
theory, one is attempting a functional integration for a theory that does not exist. One will get a result, but there 
will generically be no way of making a connection with the underlying classical theory. 

A point to consider is if one could not simply adopt in the quantum context the same attitude as was adopted in 
numerical relativity. Could one deal with inconsistent discrete theories that are just approximations of a consistent 
theory that exists in the limit? This would be a point of view that is tenable if one had a well defined way of taking 
the limit in which discretizations are infinitely refined. This is not currently available for general relativity and it is 
generally thought that such limit will be highly non-trivial given the non-renormalizability of the theory. Current 
thinking favors viewing the discrete theory as existing at a fundamental level and never taking the limit (see for 
instance [|l0| . In such a case one needs to take the issue of consistency of such a theory seriously. Even if one wished 
to consider some sort of "continuum limit" (for instance summing over all possible discretizations), it is our belief 
that one has a better chance of defining such a limit if one is obtaining it in a process where one considers a sequence 
of well defined, consistent, theories. 

In this paper we would like to discuss the canonical treatment of discrete systems. The subject is quite non-trivial. 
Up to this paper the state of the art of the subject was such that even the treatment of simple mechanical systems 
with constraints was incomplete. In particular, the literature only covers systems with constraints that are holonomic 
(function of half of the variables of the phase space). In this paper we extend the canonical formulation to include 
arbitrary constrained mechanical systems. A key new ingredient is to introduce a symplectic structure on the lattice 
based on discrete canonical transformations that implement the discrete time evolution. We then turn our attention to 
the treatment of gauge field theories on the lattice. To our knowledge, this is the first time that a complete canonical 
analysis of gauge field theories discretized on the lattice is presented, in the sense that we discretize the action of the 
theory and derive the canonical theory from the discrete action. The final results for Maxwell and Yang-Mills theories 
on the lattice coincides with the traditional transfer matrix methods, but the intermediate steps are quite non-trivial. 
We also present a canonical treatment of BF theory on the lattice. In this case, up to present no lattice treatment 
was present at all in the literature. The formulation we construct is readily applicable to the case of general relativity, 
which we will discuss in a separate publication. 

This paper is organized as follows. In section II we discuss the quantization of mechanical systems with discrete 
time evolution. In section III we extend the analysis to lattice Yang-Mills theory. In section IV we discuss discretized 
BF theory. We end with a discussion of implications for quantum gravity. 
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II. MECHANICS WITH DISCRETE TIME 



A. Summary of previous results 



Some of the difficulties of constructing Hamiltonian formulations for discrete theories can be illustrated with the 
simplest mechanical systems if we assume time is a discrete variable. The subject of "discrete mechanics" has been 
considered several times in the literature. The broad subject of "symplectic integrators" also touches upon these issues. 
However, all discussions we have found in the literature refer to systems with holonomic constraints, for reasons that 
will soon become apparent Here we will present a a discussion that includes anholonomic constraints, as the 

ones in gravity. 

Discretizing time in mechanics immediately leads to difficulties with conservation laws. Let us start by considering 
the action of a particle of unit mass in a potential, written in a discrete approximation in which time is split into 
intervals of duration e, 

N 

S = $>ife+i - 9i) + [K 2 /2 + V( qi )] e. (1) 

i=l 

If we work out the equations of motion of this action, one gets, 

Qi+i - Qi = epi, Pi - Pi-i = -V'(qi)e. (2) 

One immediately recognizes a discrete approximation to Hamilton's equation. There is an asymmetry in the equations 
in the sense that the discretization of the time derivative of q is "centered" between i and i + 1 whereas the one in 
p is centered between i — 1 and i. This theory recovers in the continuum limit the usual classical mechanics of a 
particle in a potential. However, it is well known that in the discrete theory, energy is not conserved. It is immediate 
to realize that similar problems will arise in constrained systems: the discrete evolution equations will fail to preserve 
the constraints, at least if they are anholonomic. 



A first attempt to solve the problem of non-conservation is due to T.D. Lec 12 1. His proposal consists in "pa- 
rameterizing" the system by including the time interval as a new variable. The action can therefore be written 

as, 

S = *^2p n (q n +i - q n ) - H(p n ,q n )(t n+ i - t n ) (3) 

n 

and we consider as variables pi, qi and U. When one varies with respect to U one immediately gets that H(pi,qi) = 
H(pi—i, <7i-i). This equation, with the two other equations of motion determine qi,Pi and ti 1 that is, the time interval 
is now determined by the equations of motion. 

T.D. Lee's construction merits several comments. First, notice that the use of the parameterization immediately 
suggests how to handle the problem of non preservation of constraints in general relativity. The idea is to consider 
the Einstein evolution equations together with the constraints as equations for the metric, extrinsic curvature and 
the Lagrange multipliers (the lapse and shift). One has therefore four additional equations (the constraints) and four 
additional variables. One can therefore construct a discrete system of equations that preserves the constraints in 
time. We are currently studying the feasibility of implementing this scheme in a numerical simulation of spherically 
symmetric space-times. 

There has been considerable discussion in the literature on constructing quantities that are kept invariant by the 



evolution implied by the discretized evolution equations. Maeda |13| constructed a modified Noether theorem. Logan 
fl4| proposed a method for constructing invariants not associated with symmetries of the action. Jaroszkiewicz and 
Norton Jl5| discussed the quantization of systems by constructing evolution operators. 

A whole separate chapter could be devoted to the subject of "symplectic integrators" [fl6|| . The focus of this area of 
research is usually numerical. People have noticed for a long time that it is desirable for discretization of dynamical 
equations to preserve certain conserved quantities of the classical theory. In particular the symplectic structure, or, 
as in the case of celestial dynamics, angular momentum and energy |l7fl . This is a highly developed area of research, 
in which canonical techniques have been widely used in the treatment of discrete systems. In fact, our current paper 
can be viewed as an extension of certain symplectic integration results to the realm of constrained systems, with an 
eye towards application to general relativity. 
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B. Canonical formulation and quantization of simple discrete mechanical systems 



The classical canonical treatment of discrete systems has been considered in the literature (see for instance M] ) . The 
usual treatment considers unconstrained systems and implements time evolution through canonical transformations 
of type 2 EM). The discrete evolution produces time histories of the variables. In the continuum limit, some of these 
evolutions (the ones with good continuum limit) will be solutions of the continuum equations of motion. Since we are 
interested in constrained systems, we will find that it is more convenient to work with canonical transformations of 
type 1. 

We start by denoting L(q, q) the Lagrangian of the continuum theory. We discretize time in equal intervals t n+ \ — 
t n = e and we label the generalized coordinates evaluated at time t n as q n . We define the "discretized" Lagrangian 
as, 



where 



The action can then be written as, 



L(n, n + 1) = L(q n ,q n+1 ) = eL(q, q) (4) 



q = q n and g= . (5) 



N 

S = ^L{q n ,q n+1 ) (6) 

71=0 



and its Lagrange equations of motion, 

dS dL(q n -i,q n ) dL(q n ,q n+1 ) 



dq n dq n dq n 



0. (7) 



As can be seen, in the discrete theory the Lagrangian is not a function of the manifold of g's and its tangent bundle 
q, but it is a function of q n at the n — th level and of q n +i at the next level (or cquivalcntly q n and q n ~i)- 

We now wish to introduce canonically conjugate momenta and show how the discrete evolution can be represented 
as a canonical transformation in phase space. The partial derivatives of L and L are related by, 

dL(n, n+1) dL 

dq n +\ dq 
dL(n, n+1) dL dL 
dq n dq dq 

In the continuum, the momentum canonically conjugate to q is defined as 

<"» 

and due to (jq) we define the canonically conjugate momentum in the discrete theory as 



dL(n, n+1) , , 

Pn+l = '-■ (11) 

dq n +l 



The Lagrange equations in the continuum can be written as 

dp dL 



dt dq (12) 



with p defined as in (p^|). To discretize this last expression we make the substitution 



Pn+l Pn / 10 \ 

P^ ; (13) 



and use (§), (§) and (0) to get 
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Pn+1 - Pr, 



dL _ 1 
da e 



dL dL(n,n + l) 



dL(n,n + l) 8L(n,n+l) 



-Ln+1 



Pn+1 + 



dL(n, n + 1) 



Then 



dL(n, n + 1) 
dq n 



Summarizing, the discrete Lagrange equations are, 

dL(q n ,q n+1 ) 



Pn+1 = 



dq 



Pn 



n+1 



dL(q n ,q n+1 ) 
dq n 



(14) 



(15) 



And it should be noted that these equations define a type 1 canonical transformation from the variables (q n ,Pn) to 

(q n +i,p n +i)- 

The variables (q n ,Pn) constitute a phase space and as usual, if the system of equations ( |l5| ) is non-singular one 
can introduce new coordinates (g n +i,p n +i) for the phase space. This transformation is canonical, in the sense that it 
preserves the symplectic structure, 



{q n ,Pn} = {q n +i-,Pn+i} = 1. 



(16) 



The generating function (we use the terminology o f |l9| ) for the transformation is F(q n ,p n ) — —L(q n , q7i+i(q n ,Pn)), 
where q n +i(q n ,Pn) is given by inverting the first of (|15|); it is immediate to check that it satisfies the equations of a 
generating function, 



OF 

dq n 
dF 

dp n 



ln+1 



-Pn+1 



dq n ^ 



dp r , 



-Pn+1- 



(17) 
(18) 



One can now introduce the generating functions of canonical transformations of the various types that correspond to 
discrete time evolution. If the variables (q n , q n +i) are used for coordinates of the phase space, then one can introduce 
a generating function of type 1, 



Fi(q n ,q n +i) = -L(q n ,q„+i) 

Pn = ^ = {^l,Pn} 
oq n 



Pn+1 = 



8Fx 



ln+1 



= - {Fl,Pn+l} ■ 



(19) 
(20) 

(21) 



It is then possible to define generating functions of type 2,3 and 4. For instance, if the phase space is coordinatized 
by {qn,Pn+i), solving the first of ( |l5|) for q n +i one can introduce a generating function of type 2, 



F 2 (qn,Pn+i) = q n +iPn+i - L(q n ,q n+X ), 
dF 2 



Pn 

qn+i 



{F 2 ,Pn} : 



dF 2 



{F 2 ,q n +i} 



(22) 
(23) 

(24) 



dpn+l 

To facilitate the comparison with the continuum limit case, we can introduce a function H(p n+ i, q n ) given by, 

F 2 (q n ,p n+ i) =p n+X q n + H(p n+1 ,q n ), (25) 
such that one recovers the discrete Hamilton equations, 

dH 



qn+i = qn 



Pn = Pn+1 



dpn+l 

dH 

dq n 



(26) 
(27) 
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The quantization of these finite canonical systems is achieved by representing the finite canonical transformations 
that materialize the time evolution as unitary operators. We illustrate this with a simple example, a particle in a 
potential. The Lagrangian is given by 

L{q n , q n +\) = m ^" +1 2e V(q n )e (28) 

the canonical momentum is given by p n +i — m(Qn+~L ~ <ln)/ e i from which we can get q n +i = Pn+x e /' m + <?«.• The 
generating function is given by, 

P 2 

F 2 = Pn+lQn + -7T^- e + Vfan+i)* = Pn+iQn + H(p n+1 ,q n ), (29) 
2m 



and the equations of motion (taking into account (^3 24)) are 



i P n + 1 /on\ 

q n +i = q n H e (30) 

m 

Pn = Pn+i + V'(q n )e (31) 



which can be solved for p n +\ as, 



q n +\ = q n + — e - V'(q n )— (32) 
m m 

Pn+i =Pn - V'(q n )e. (33) 

We now proceed to quantize the system. We choose a polarization such that the wavefunctions are functions of 
the configuration variables, ^/(q n +i)- The canonical operators have the usual form. The evolution of the system is 
implemented via a unitary transformation. An immediate computation shows that the operator that implements the 
transformation p n +i — Up n U\ q n +i = Uq n U> is given by, 

^ = exp ^__jexp(^— j. (34) 

At a quantum mechanical level the energy H°(q n+ i,p n+ i) is not conserved, as we expected from the fact that it 
was not conserved classically. It is remarkable however, that one can construct an "energy" (both at a quantum 
mechanical and classical level) that is conserved by the discrete evolution. To construct this energy we use the 
Baker-Campbell-Hausdorff formula, 

exp (X) exp (Y) = exp (x + Y + ± [X, Y] + -L ([X, [A, Y]] + [Y, [Y, X}}) +..?), (35) 

and one can therefore write U = exp (^j-H cS (q i ,p i )') where H c g = H°(q,p) + O (e 2 ), which is immediately conserved 
under evolution. It is straightforward to write down a classical counterpart of this expression. 

C. Canonical formulation and quantization of constrained discrete mechanical systems 

We now extend the analysis of the previous subsection to the case of a constrained system. Here it will become 
clear the advantage of using canonical transformations of type 1. We start with an action written in first order form, 

L(n,n+ 1) =p n (q n +i - q n ) - eH(q n ,p n ) - \ nB 4> B {q n ,Pn) (36) 

where we assume we have M constraints B = 1 . . . M. 

We now exhibit the construction of a type 1 canonical transformation. We construct the appropriate canonically 
conjugate momenta using the first of (fl5|), 

_ dL(n,n + l) _ 

r n+l — a — Pn {O I ) 

oq n +i 

p p _ dL(n,n+l) _ ( 
P " +1 " d Pn+1 ° (38) 



)A b _ dL( n,n+l ) 



p^i= a \ =a (39) 



G 



To determine the equations of motion for the system we start from the second set of equations (|l5|), 

dL(n,n + l) dH(q n ,p n ) dcf) B (q n ,p n ) 

Pn = 5 = Pn + e ^ h * nB 5 ( 40 

oq n dq n dq n 

dL(n,n + l) dH(q n ,p n ) d(j) B (q n ,p n ) 

r n — q — -\q n +l - In) + e o r A nB T. l 41 J 

Opn Cfpn Opa 

Pn" =<t> B {ln,Pn)- (42) 
Combining the last two sets of equations we get the equations of motion for the system, 

dH(q n ,p n ) d(j) B (q n: p n ) 

Pn-Pn-i = -e 5 Kb 5 (43) 

oq n aq n 

dH(q n ,p n ) d(f) B (q n ,p n ) 
q n +i-q n = £ 5 h A„ B 5 (44) 

S (9n,Pn)=O. (45) 

Superficially, these equations appear entirely equivalent to the continuum ones. However, they hide the fact that in 
order for the constraints to be preserved, the Lagrange multipliers get fixed. Another way to see it, is that = p n 

and therefore it is immediate that the Poisson bracket of the constraints evaluated at n and at n + 1 is non- vanishing. 
We then consider the constraint equation (|45|) and substitute p n by (|37j), 

B (q n ,P* +1 ) = O. (46) 

We then solve (Q) for q n and substitute it in the previous equation, one gets, 

</> B (g„ +1 ,P n 9 +1 ,A„ B ) = 0, (47) 

and this constitutes a system of equations. Generically, these will determine 

KB=K B {qn+uP q n+1 ,v a ) (48) 

where the v a are a set of free parameters that may arise if the system of equations is undetermined. The eventual 
presence of these parameters will signify that the resulting theory still has a genuine gauge freedom represented by 
freely specifiable Lagrange multipliers. 

The final set of evolution equations for the system is therefore given by (^ 44 ) where the Lagrange multipliers are 
substituted using ([48|). 

At this point it is worthwhile discussing the physical meaning of the determination of the value of the Lagrange 
multipliers. This is an issue that is potentially confusing. Lagrange multipliers are normally associated with gauge 
symmetries. The fixation of the value of the Lagrange multipliers seems to suggest that the constructive procedure of 
the discrete theory somehow has selected a preferred gauge for us. A priori this might appear very unnatural. How 
could the procedure know which gauge to choose? A point to emphasize is that in this section we have absorbed in the 
definition of the Lagrange multipliers a factor of e (the evolution parameter interval) when we replaced the integral 
of the action by a discrete sum. What gets determined is therefore "A x e" . For a completely parameterized theory, 
since there is no explicit reference to e, one can choose it arbitrarily and therefore redefine the value of the multiplier 
arbitrarily. A different way of putting this, which might help later understand the situation in general relativity, is to 
consider that once the lattice has been established, the gauge is fixed, but one has still the freedom of where to place 
the lattice in an arbitrary manner. 



III. YANG-MILLS AND MAXWELL THEORIES ON A DISCRETE SPACE-TIME 



Having set up the canonical quantization of mechanical systems discretized in time, we now turn our attention to 
field theories. In this case we also wish to discretize space. Although the words "Hamiltonian formulation" are used 
frequently in the lattice context, the formulation is not constructed in the traditional canonical sense. A popular way 
of obtaining a Hamiltonian operator for quantized lattice gauge theory is the "transfer matrix" method. This consists 
in discretizing the space-time action and reading off from the path integral the operator that would correspond to the 
infinitesimal evolution operator in the limit in which the spacing in time goes to zero. There is no construction of a 
classical discrete theory that is later quantized. Even the definition of canonical variables for such a classical theory 
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is problematic and has received some attention |20[|. What we will do in this paper is precisely that: we will construct 
a classical canonical discrete lattice theory for Yang-Mills, that we will quantize and we will show that we recover the 
same results as those of the usual transfer matrix approach. 

We will first explore the construction of Yang-Mills lattice theory with a continuum time parameter. Although we 
know that ultimately this is not what we want, since as we argued it is unlikely a formulation like this will be useful in 
the gravitational case, it will prove instructive for later constructing a theory discretized in space and time. We will 
then consider the quantization of Maxwell theory on a discrete space-time. Maxwell's theory has all the ingredients 
of interest and is simpler than Yang-Mills theory. 



A. Discrete lattice Yang— Mills with continuum time 



We consider a hypercubic lattice in four dimensions with spatial spacing a and spacing in time ao- We label the 
vertices of the hypercube with an integer i. Each link connecting two sites is identified as a pair {i,j}- Associated 
with each (oriented) link there is a holonomy Uij, we shall assume for simplicity that it takes value in SU(2). The 



holonomy of the link traversed in the opposite direction is given by the inverse Uji 
by, 



S 



aoa 



E 



5 P 
.4 



2 Tr((7p) 



U ij . The Wilson action is given 



(49) 



where the sum is over all elementary plaquettes P and Ap is the area of the plaquettes. Sp is —1 if the plaquette is 
spatial and +1 if it is timclike. Notice that we are considering the Lorentzian action, our treatment does not require 
a Wick rotation to the Euclidean action. We now will take the limit in which the spacing in the time direction, ao 
goes to zero. In order to do this we will relabel the links with two indices I, n where the index n labels the position 
of the vertex "in time" and the index I runs through all the oriented links that form the spatial sub-lattice L n at a 
given time n. We denote by »l and I* the vertices at the origin and end of the link I. We rewrite the action and we 
make explicit the elementary plaquettes that are completely spatial and those that have temporal links, 



S = - 



E 



i 



Pn 



where holonomies along timelike links are denoted by the spatial label of the vertex they start and end at, and the 
time level they start at, like Vi m<n . The second term in the Lagrangian is a sum over all spatial elementary plaquettes 
at time n, P n . We introduced the following notation for the components of the holonomy, XJ\ = U^T 1 where Uf are 
real, the T are given by T° = I/y/2, T a = -ia a /V2, a = 1 ... 3 and a are the Pauli matrices. 

We will now take the limit in which ao — > 0. The holonomies along the timelike direction become, 



V v , 



1 + X v ao + 



(A.) 



K^ = l-A„a + (A„) 2 ^ + 
a o " 

Ul n +\ = Uln + aoUln + — Ul^ n 



(50) 
(51) 
(52) 



where v is any vertex of L n and A„ = \%T a with A° real. We will use as generalized coordinates the variables Uf and 
\ T V , normalizing the Lagrangian in such a way that S = aoL(n, n + 1) 
The Lagrangian therefore becomes, 



L 



a E Tr 

\l\EL n 



1 



A; 



\2 



+^E^) 



3 

a 3 V ai ^ 



i 



1=0 



(53) 



We have dropped the subscript n since all variables are now evaluated on the same time slice. The above sum goes 
over all links traversed only once in a given direction and Uj = U^ 1 . We have also included an additional term in the 
Lagrangian, which is a constraint that fixes the holonomies to take values in SU(2) provided the A^,'s are real. 
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It can be seen that the discrete Lagrange equations of motion coming from this Lagrangian lead, in the limit 
a — > to Yang-Mills equations. The discretization in this case does not break any gauge symmetry and therefore the 
discretized theory is automatically consistent. 

We now proceed to work out the canonically conjugate momenta, 



u ui t i 



n? 



dL t q dL 
dU? 8U, 
dL 
dai 
dL 



= 



Performing the Legendre transform, the Hamiltonian becomes, 

' ' ' ' ^E^ 



H= £ Tr 

\i\eL n 



2a 



nfnf + nY(u l x l .-x.iU l ) 



« 3 E«'E[(^) 2 -2]+ E wnr+EE^W 



1=0 



\l\GL„ 



(54) 

(55) 
(56) 



(57) 



where p\ and fx va are Lagrange multipliers. We now need to ensure that the evolution preserves the constraints. 
Ensuring the conservation of (pq) leads to 



E N 



- 2 



7=0 



= 



(58) 



and conservation of (|56j) leads to Gauss' law, as we will soon discuss, 
and of Gauss' law. Conservation of (|5^) leads to 



We now need to ensure the conservation of <\5 



Tr 



0. 



(59) 



This implies that the product E\ = UiHf is an element of the su(2) algebra and is purely imaginary. To understand 
this, we need to notice two properties. First of all, the object UiU-Y transforms under gauge transformations as a 
gauge vector at the point •/. The second property is that Gauss' law can be rewritten as, 



\i\eL n 







(60) 



where k v is a smearing function. This is the usual expression of Gauss' law in lattice gauge theory, and shows that 
the quantity UIl u plays the role of an electric field, and the Gauss law just states that the outgoing flux of electric 
field lines from point i vanishes. Conservation of Gauss' law is guaranteed, since the action is gauge invariant, but it 
can be explicitly verified. 

We still need to conserve in time the constraint that states that UH U is a purely imaginary element of the algebra. 
The resulting equation determines the Lagrange multiplier a;. The fact that the Lagrange multiplier is zero implies 
that there will be second class constraints. Indeed, constraints (|5|,|5^) are second class. 

To handle the second class constraints we can introduce Dirac brackets. The Dirac brackets become, 



{ml . iuvfc) - (#*g- {vl) AC m BD ) 



(61) 



This relationship was derived, via an independent constructive method by Renteln and Smolin p0| |. Based on this 
Dirac bracket one can work out the Dirac brackets of the electric field with the holonomy and with itself and one 
recovers the usual brackets that people assume in lattice gauge theory and that approximate in the limit a — > the 
Poisson bracket in the continuum of the electric field with a holonomy. 

We now proceed to write the Dirac extended Hamiltonian, which is achieved by imposing strongly the second class 
constraints and adding the first class secondary constraints we found. Performing the sum, we get, 
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H=— {Ef + «.,£?, + Kl .Ej) ~ — £ Tr (U P ) . (62) 

|{|ei„ f 

And we see we recover the well known Hamiltonian of Yang-Mills theory on a lattice. This example is revealing in 
that if we had simply taken equation ( |53| ) and attempted a direct quantization ignoring the fact that there were second 
class constraints it would have been very difficult to reproduce the usual Yang-Mills Hamiltonian in the lattice. All 
attempts to quantize the Hamiltonian constraint of general relativity have ignored the complex canonical structure 
of the discrete Hamiltonian being quantized. To address this issue is the main point of our approach. 

This example also exhibits how treating time as a continuum variable while discretizing in space yields a complicated 
canonical structure for the resulting theory. In the case of general relativity it would be expected that complications 
will be even greater, since the constraint algebra intimately ties space and time. We will therefore from now on only 
discuss field theories in which both space and time are discretized. We start by considering the example of Maxwell 
theory. 



B. Maxwell theory in a discrete space-time 

As in the case of Yang-Mills we start by writing the Wilson action, particularized to the case of a U (1) connection, 
L(n,n+ 1) = ~— ]T {Ul n+1 V: hn U l>n V lt>n + U^+.V.^U^V^) + ]T U Pn 

|I|6i„ Pn 

+a a 3 J2 ( U i^ U tn - 1) + a o« 3 Pv.n (V v , n V* n - l) (63) 

\l\EL„ v£L n 

where as before L n es the spatial sub-lattice at a given time n and we label the links with I, the vertices with v and 
the plaquettes with p n . The sum |Z| € L n goes over all links traversed only once in a given direction. The holonomies 
along spatial links are denoted by U^ n where / is the link and n the "time level" at which they live. Holonomies along 
timelike links are denoted by V v>n where v is the spatial label of the vertex they start in. In this case U and V are 
complex numbers. Given a complex number A we will denote A and A 1 its real and imaginary parts respectively and 
we will write A = A I T I where T° = 1 y T 1 = i. We will take as generalized coordinates U 1 , V 1 and the Lagrange 
multipliers a^ n y 

We shall now derive the canonically conjugate momenta and the equations of motion for each variable. For a and 

/3 we get 

dL(n,n+l) 

+1 = — = (64) 

oai <n+1 

dL(n,n+l) 3 - , 

n ',» = a = _a ° a K U i^ U i,n - !) 65) 

< n+1 = d J£p±ll = (66) 
Kn = - dL ^ H+1) = -a a 3 (V v , n V* n - 1) . (67) 

Equations (^) y (|6^) are primary constraints and their preservation leads to the secondary constraints 

Ui,nUl n -1 = (68) 
V v , n V* n -1 = 0. (69) 

And these constraints are automatically preserved. 

The momentum canonically conjugate to the variable J7/, n +x is 

n^ +1 = ^?;; +1) (^r (70) 

OU l,n+l 



and from the Lagrangian (]63l) we get 
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u 

l,n+l 



-v, 



aa 



PL 



(71) 



where Vpi_ n = t^.nVj.^t/^^V^ n is the holonomy along the plaquette that starts at the link I and continues along 
a timelike link. It should be noted that the evolution equations together with the constraints that imposed unitarity 
on the holonomics imply that the momenta are unitary up to a factor, 



L1 l,n 



(72) 



which we will later need to impose strongly as a constraint. 
For the momentum 11^ the Lagrangian leads to 



UirJIi _ = -2a a ai, n H Vpi. n U Pl 

u Pi 



(73) 



where Up l is the holonomy along a spatial plaquette containing the oriented link I and the sum is over all these 
plaquettes. The real part of this equation determines the value of the variable a^ n and the imaginary part leads to 
the following evolution equation (on the constraint surface) 



Ul,nW 



(74) 



For the variable V the canonically conjugate momenta at times n and n + 1 are 



Kn+l = 



V IF 



-2a Q a (3 V , 



^ (Sv.lmVpi 



\i\eL n 



(75) 
(76) 



Studying the evolution of the constraint ( f75| ) , the real part of equation ( |76| ) determines the value of the variable 
and the imaginary part (substituting the previous constraints) leads to the constraint that will end up being 
Gauss' law, 



E 

\i\eL T . 



CLv 







(77) 



where ai lV = 5 V: ,i — S Vi i». Notice also that there is no evolution equation for the V J s, we will later see that they are 
free parameters in the theory that correspond to gauge transformations. 

Let us review the constraint structure of the theory. We have the following constraints, 

(78) 
(79) 

(80) 

(81) 
(82) 
(83) 
(84) 

\i\eL n 





Ul,nUi n 


= 1, 




V V* 


= 1, 






_ a 2 
a l 


v v 


ttV . 

n+l iJ -u,n+l ' c - c - 


= 0, 


v v 


yrV 

n+l 1J -D,ri+l c ' c - 


-o, 




ll l,n+l 


= o, 




11 v,n+l 


= 0, 


[Ul,n+lKY,n+l ~ c - c ] 


= 0, 



into the two linear combinations (81,8^). The set of constraints 
I) is first class. The remaining constraints are second class. The constraints that state that U and II are unitary 



where we have divided the constraint n^„ +1 



are second class. We need to work out the non-vanishing Dirac brackets, given by 
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( n ^ 



R-i€} Dlrac = 2 - , rr v ':f l " y TI r (86) 



{^,n,W} D . - - 2 W^ (87) 

L JDirac (n[>,„) + nff n c/ i ,„ 



{K,n,n^„} Dirac - 1 

\v v , n ,(nl n )*} =-(K,n) 2 (89) 

I* J Dirac 

and their complex conjugate. At this point we are, in principle, finished with the classical treatment. What we now 
need to do is to perform a quantization, which entails finding the physical space of states (states that are annihilated 
by the first class constraints) and a unitary evolution operator that embodies the finite canonical transformation that 
yielded the above equations of motion. 

In order to identify the evolution operator it is desirable to simplify the structure of the theory. This can be achieved 
through the introduction of real variables 

^^(^-(n^^r), (90) 
(n^T4, n -(n^)*(K,„)*), (9i) 



F v = - 

v ' n ~ 2 



that have simple Dirac brackets, 



{Ei, n ,Ui, n } Biiac = -%Ui, n (92) 

n£„} Dirac =+*n,% (93) 

«„,V , n } Dirac =-iV 0lB (94) 

{<„,n£ n } Dirac =+in£„. (95) 



In terms of these variables, the first class constraints ( p2| ) and ( |85| ) become, 

<™ = 0, (96) 
J2 <ri.vEi,n = 0, (97) 

\l\GL n 

and we recognize in the last expression Gauss' law. The simple canonical structure of these variables immediately 
suggests the introduction of variables canonically conjugate to Ei t7l and E^ n , which we call ipi tn and <p v ^ n respectively 
and are given given by, 

U l>n = exp (iipi t n) (98) 
V v>n = exp (i<fi v<n ) (99) 

with canonical Dirac brackets 

{ V i,„,E lin } Biiac =1, (100) 
Kn,<„} Dlrac =1- (101) 
The new variable has remarkably simple evolution equations. Starting from (|74|) we get 



Ei, n+1 = E lin + i^Y] [U(Pt) - U*(P t )] . (102) 
We need to find an evolution equation for ip. In order to do this, we consider the definition of E, equation (|9 



evaluated at n + 1, and we substitute nP n+1 using equation (71). The result is an alternative expression for E^ n+1 , 
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2a 

We now rewrite the right hand side in terms of the angular variables, 

£7.71+1 = Sin ((fl, n +l - 4>lm,n ~ fl,n + <j>»l,n) 

a 



If we now equate the right hand sides of (102) and (104) we get the desired evolution equation, 

= <f>l;n + <Pl,l 



(103) 



(104) 



(105) 



With the construction of the simplified variables, we are ready to quantize the theory. We have a canonical 
Dirac bracket among variables, so it is immediate to pick a representation for operators and wavefunctions. We 
choose wavefunctions ^{tpi, (p v ,ai, f3 v ). Imposing the first-class constraints (|82]-85) implies that the wavefunctions are 
independent of a, /3, 0, and that as functions of ip they have to be gauge-invariant (this is the content of the quantum 
Gauss law). 



The quantum evolution equations are obtained by promoting all quantities in equations (|102| ) and ( 105 ) to quantum 
operators. The important task is to construct a unitary evolution operator that would embody the time evolution 
implied by the above discrete quantum evolution equations. We proceed in steps. In order to reproduce the evolution 
equation for E we notice that, 



E hn+ i = exp I ^2 [Upi + Up} J Ei, n exp I ^ \u P i + U P l 
\ a pi ) \ a pi 



AEi^A- 1 



where by 



A = exp —i 



(106) 



(107) 



we denoted the "partial" evolution operator y we denote Uj = U, 1 . To reproduce the equation for ip we notice that, 



»./,n + sax — Ei i7l+ x 



v a 



(108) 



and therefore, 



I'm.n ) E, 



— 1 / a - 

(pin + sin — Ei. 

V a 



<pi, n +i = A exp ( i ^ (<f>'l',n 

\l'\£L„ 

i E {$*V,n _ 4>l';nj Ey m J A^ 1 

\l'\&L n J 

A exp I i ^2 (4>»i',n - </V.,ri) Ei^ n + J dxsin -1 



x exp — % 



|i'|eL„ 
x exp I —i 



We have therefore obtained the evolution operator, 
.00 



i ^2 (J>»l',n - 4>V;nJ El' ,n + I dxSXLX 1 ( — X 



x=B,i 



x=E,, 



x exp li 



"f>ti', n Ei', n + ^ I dxsin 1 (—x 



<Pl,n 



(109) 



(110) 



Where the sums are over all possible plaquette and link orientations and we are using that Ej — —Ei. This evolution 
operator coincides with the one obtained with the "transfer matrix method" in ordinary lattice Maxwell theory (the 
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term with the Gauss law is not recovered since the transfer matrix method is usually applied in a gauge fixed context). 
See for instance formula (5.13) of reference [pl| . One can take the continuum time limit of this expression ao — > 
and one recovers the usual Hamiltonian of lattice Maxwell theory. 

This concludes our discussion of ordinary lattice gauge theory. It is evident that the proper use of discrete canonical 
techniques is possible and reproduces well known results of the transfer-matrix method. This method is geared towards 
theories with a definite time variable and an associated true Hamiltonian. We will now discuss fully constrained 
theories, which are less suited for treatment with transfer-matrix methods. 



IV. BF THEORIES ON THE LATTICE 



To our knowledge, the transfer-function method has not been applied to fully constrained field theories (like BF 
theory or general relativity) without performing a gauge fixing (for gauge fixed attempts to lattice gravity see p2|). 
The techniques we are introducing in this paper, however, are not limited in their ability to handle fully constrained 
theories. This is good, since it gives hope that they could work in the gravitational case. To begin with we will consider 
the case of BF theory. BF theory has been quantized by many techniques, since it has a finite number of degrees of 
freedom. The interest of our approach is that it will yield a discrete theory that has a proper canonical formulation 
and that contains within it solutions that correctly approximate those of the continuum theory at a quantum level. 

The treatment we discuss applies to BF theories in any dimension. We consider a continuum action of the form, 

S = J Tr (£> A F) (111) 

where F is a 2-form associated with a connection (to make things concrete, let us assume it is an su(2) connection, 
although our method extends easily to other semisimple groups) and B is a D — 2 su(2)-valued form. 

In the continuum, if we perform an N + 1 decomposition, where N = D — 1, one ends up with a system of first-class 
constraints given by, 

D a E a = 0, (112) 
F ab = 0, (113) 

where E a — e abc 'Bb c ... is the spatial dual of the spatial pull-back of the form B and F a b is the spatial pull-back of the 
curvature F, indices a,b, . . . run from 1 to D — 1. So the content of the theory is that the connection is flat and there 
is a Gauss law. Quantum mechanically, the constraints are solved by functionals of connections that have support on 
flat connections only and are SU(2) invariant. The theory only therefore has a finite number of degrees of freedom 
and is non-trivial only in non-trivial spatial topologies. In (2 + 1) dimensions the theory (with an SU(2) group) is 
equivalent to Euclidean Einstein gravity. 

We discretize the BF action in the following way, 



L(n,n+1) = ^Tr 



B n , v h n v h n vJrei {h n v+e2 ^ {h n ,v) + ^n,v^n,vh n+ i_ v {Vn^+ex ) (^n,ti) (H4) 



z 

B n , v Vn,vh n+ i v (Vn,u+e 2 )^ i} l n,v) ""(" ^n.v {V n _ v V^ v — /) + ^ ] ^n.v f^n,D (^n,u) — 



k=l 



To explain the notation we refer to figure IV . From now on we will assume we are in 2 + 1 dimensions in order to 
simplify notation although it will be evident that simple generalizations of the formulae will hold in any number of 
dimensions. By h\ v we mean the holonomy along the direction k (in the 2 + 1 dimensional case it is either along 
the elementary unit vectors e\ or e-i) starting at the lattice point labeled by the time step n and the spatial point 
v (in 2 + 1 dimensions v is labeled by a pair of indices. By the index v + we mean the lattice point arrived to 
by following starting at lattice point v. The B fields live D — 2 hypersurfaces dual to the plaquette on which 
we compute the holonomy representing the F and are elements of the algebra of SU(2). In 2 + 1 dimensions B is 
a one-form we therefore write it with one index. The vertical holonomies V n<v are only labeled by the lattice point 
they start at. We will assume that the holonomies are matrices of the form h — ^T 1 , V = V I T I where T° = I j\[2 
and T a = —ia a / \[2 where a a are the Pauli matrices. From now on when we take variations we will assume that one 
has performed the variations with respect to the components h 1 and reconstituted the remaining equations as matrix 
equations. A and fi are Lagrange multipliers that enforce, given the above condition of the holonomies, that they are 
elements of SU{2). 
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(n,v) 



1 (n,v+e ) 



FIG. 1. The elementary cell which we use to discretize BF theory. 

The fundamental canonical variables of the theory are therefore /j„'^, V£ v , -B„S) v and /J, n>v . The equations of 
motion for the variables v are, 



n 



n 



a" 



dL(n, n + 1) 

UA n+l,v 

dL(n, n + 1) 
d A£ ., 







= -Tr 



(115) 
(116) 



Preservation in time of the primary constraint (|115|) leads to the secondary constraint 



Tr 



h k (h k ) f - I 

n n,v \"'n,vl 1 







(117) 



and preserving in time this constraint does not yield any new constraints. The equations of motion for the variables 
f^n,v are 



TP — n 

11 n+l,v — u 

LP = -Tr 



v niV (v nv y - 1 



Preserving in time constraints(118) leads to the constraints 



Tr 



V n , v {V n J - I 



= 



(118) 
(119) 

(120) 



which is preserved in time without introducing any new constraints. 
The equations of motion for the variables are 



n 



dL(n, n + 1) 



n+l,v 



dB 



k . a 



U B = - 

n.v 



(T a y = o 

n-\-l,v 

dL{n,n+l) f 



dBnv 



(T a Y = -Tr [htf T a ] (T a ) 



(121) 
(122) 



where {i,j, k} is a cyclic permutation of {0, 1,2} and h l n \ is th e ho lonomy starting from n, v a round th e elementary 
plaqucttc in the plane i — j. Preserving in time the constraint (121) together with constraints (117,120) leads to 



hi 3 — J-3 I 

n,D w n,v 

with a l n J v — ±1. This sign is arbitrary and can change from plaquette to plaquette. 
The equations of motion for the variables V„ v are 



(123) 



v _ dL(n,n+l) , A t 

LL n+l,v — |J 







V T\ v — /i u 1 R 
v n,v LL n : v ~ ll 'n,v £} n,v 



01 D 2 .02 nl —(f)' 1 \^ h 10 R 2 h 1 

' n,v n,v \ n,v—eij n,v—ei n,v—ei n,v—ei 

( h 2 \t .20 R l 1,2 , 9 r 

\ n n,v — e 2 l "'n,v — e 2 JD n,v — e 2 "'n,v — e 2 ' 1 ■ 



(124) 



(125) 
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Preserving in time the constraints (124), together with ( |123| ) leads to determining the multipliers 
and the equation of motion, 



lOI d2 _ l02 nl _ /.l \^ h 10 B 2 h 1 + (h 2 V h 20 B 1 h 2 

n,v n,v n,v n,v \ n,v — e\) ' n,v— e\ n,v— ex n,v — ei ' V n,v— e2/ n,w-ea n,u-ea n. 



The equation of motion for n^ +1 v is 



Ll tt/i 1 _ T/f /,01 R 2 T/ 

' t n+l,'U 11 n+l,t) — v n, u"n, ti- ™, « v n,v 



and using (122) we obtain the primary constraint 



Ti- 



ft 1 M hl 

"'n+l,v LL n+l,v 



(126) 



(127) 



(128) 



(129) 



Preserving in time this constraint, together with the equation of motion for 11^ v leads to the fixing the multiplier 



,, — 



and the equation of motion 



l hl 



h 1 IL h 



h 10 B 2 — h 12 B° 



111 \ t , 2 1 nO .2 

V n,f — e2 / n,v — e2 ri,v — e2 n,v — e2 



We proceed in a similar fashion for the variable h% v . The equations of motion lead to 



Tr 



h 2 n h2 — _T/t h 02 r 1 v 

""n+l,v 11 n+l,v — v n,v"'n,v JD n,v v n,v 

h 2 M h2 



K .. = 



h 2 M h2 —h 21 B° - h 20 B 1 — (h 1 ^ h 12 B° h 1 
We now consider the variable, 

E n ,v = hn,vH-n,v k € {1, 2} , 



(130) 
(131) 

(132) 
(133) 
(134) 
(135) 

(136) 



that our experience with Yang- Mill s sh ows should play the role of an electric field along the direction eu at the point 
{n, v}. Due to the constraints (129) y ( |133| ), the variable in question is an element of the algebra, i.e., is a traceless 
matrix, 



Tr = 0. 

We also define the electric field in the direction — ej, at the point {n, v} by 

E nv = h nv Yi nv — (h n „_ e . ) e ) — — (h nv _ e ) E n v _„, h" 



e k '"n,v — e k 



Using the equations of motion and constraint it is straightforward to get 

77U _i_T/t /jOl 

^n+l.u — ^ v n.v ll n,v D n,v v n,v 

7^2 Trf ^0 2 nl T/ 

-^71+1,1; — V n,v"'n,v- D n,v v n,v 

P T _ T/ f ( h l ^h Q1 B 2 h 1 V 

& n+l,t) v n,v V l n,v-ei I n n,v — ei ^n^v-ei n n,v — ei v n,v 

p2 . _i_T/t (1,2 ^t/,02 R l 1 2 v 

& n+l,-u ' v n,v\ n n,v-e 2 ) ""n,v—e 2 n,,v — e2 n,v — e2 n < v 



and from here one can show that (127) is equivalent to 



(137) 

(138) 

(139) 
(140) 
(141) 
(142) 

(143) 
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This is the familiar expression for Gauss' law on the lattice. 

We have noted that the holonomy along the plaquette 1 — 2 is proportional to the identity. Therefore one is dealing 
with a flat connection on the plaquette. As expected, the discrete theory admits more solutions than the continuum 
one, in the sense that a possible solution would be a connection that makes the holonomy +1 on certain plaquettes 
and —1 on others. Such solutions will not admit a continuum counterpart. 

We need to check that the Dirac brackets of the elementary variables lead to the expected results. Here we can j ust 
refer to the discussion we did in the Yang-Mills case, since the structure of the second class constraints (117 137| ) is 
similar to the one we encountered in that case. 

The quantization of the theory is now immediate. The Dirac bracket between the electric field and the holonomy 
is similar to the one in the Maxwell case (B2l), 



{(E*, V ) A ,h* iV } =iT A h* >v (144) 

where we have introduced the component notation in the algebra and the T A, s are the generators of su(2) we 
introduced before. 

One can now consider wavefunctions that are functions of the holonomies on the links of the lattice, ^[h^ lv ]. 
Following the Dirac method we need to impose the constraints as operator equations. Gauss' law will just state that 
the wavefunctions are gauge invariant. In this case it will imply that the wavefunctions are linear combinations of 
traces of holonomies along closed loops. 



We now need to impose the spatial projection of the constraint (123) as an operatorial equation, 

< 2 „ = <U (145) 

In the space of functions of holonomies the equation becomes multiplicative and the solution is that the holonomies 
along all elementary plaquettes should be ±1 times the identity. We therefore see that we recover as solutions of the 
theory the space of gauge invariant functions of a flat connection. Except that in the discrete theory it could happen 
that in some plaquettes the holonomies are +1 times the identity whereas in others the holonomies have a negative 
sign. Such solutions are possible in the discrete theory and do not have a counterpart in the continuum. We see that 
the discrete theory therefore has considerably richer solution structure than the continuum theory. If one however 
imposes that the solutions be continuous from one plaquette to another, one is only left with the traditional solution. 

The unitary operator that implements the finite canonical transformation is, for this theory, given by an exponential 
of a combination of the first class constraints with arbitrary coefficients. 



V. CONCLUSIONS 



We have generalized classical mechanics to systems in which space and time are discrete. Time evolution is 
implemented through a finite canonical transformation. Constrained systems can be incorporated in the scheme. We 
discussed in detail how to handle Yang-Mills theory, Maxwell theory and BF theory and encountered no problem 
setting up the classical and quantum theory. In general, the discrete theories contain many solutions that do not have 
a good continuum limit, as we have seen in the BF example, and the correct solution has to be chosen by requesting 
that it have the appropriate semiclassical limit. 

The framework we have set up is completely ready to be applied to the case of general relativity, which will have 
many similarities with the BF case. An attractive advantage is that the formalism does not require an Euclidean 
rotation and operates directly in terms of the Lorentzian signature theory. We will discuss the case of general relativity 
in a forthcoming publication. Having a sound theoretical foundation for dealing with discretized gravity will allow 
quantize the theory in a systematic way and will enable us to make connections between the canonical and path 
integral approaches. 
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